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A new paradigm for active galactic jet kinematics has emerged through detailed investigations of BL Lac objects using very 
long baseline radio interferometry. In this new scheme, most, if not all, jet components appear to remain stationary with 
respect to the core but show significant non-radial motions. This paper presents results from our kinematic investigation of 
the jets of a statistically complete sample of radio-loud flat-spectrum active galaxies, focusing on the comparison between 
the jet kinematic properties of BL Lacs and flat-spectrum radio-quasars. It is shown that there is a statistically significant 
difference between the kinematics of the two AGN classes, with BL Lacs showing more bent jets, that are wider and show 
slower movement along the jet axis, compared to flat-spectrum radio-quasars. This is interpreted as evidence for helically 
structured jets. 
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1 Introduction 

Black holes have been at the center of astronomical research 
for many years, being one of the most prominent predictions 
of the theory of General Relativity. There is a large amount 
of indirect evidence supporting that most, if not all, galax- 
ies host a supermassive black hole (SMBH) at their cen- 



ters (e.g., Blandford 1986). It is widely believed that this 
SMBH is the main constituent of activity in active galactic 
nuclei (AGN) (for a review see e.g., Begel man et al.| 1984 1. 
Given their very luminous nature, AGN are therefore some 
of the best probes that we have today to study the proper- 
ties and evolution of SMBH out to the earliest times of the 
Universe. It is believed that AGN are made up of a number 
of building blocks including a central SMBH, an accretion 
disk, an obscuring molecular screen (usually referred to as a 
torus), ionization regions producing broad and narrow emis- 
sion lines, and in some cases powerful collimated outflows, 



known as AGN jets (e.g., Antonucci 1993 Urry & Padovani| 
[I995l >. 

Jets are found in a plethora of astrophysical environ- 
ments, among which, young stellar objects, massive X-ray 
binaries, possibly pulsars, y-ray bursts, and AGN. They 
have been observed for the first time in the optical in the 
nearest active galaxy to us, M87 ( |Curtis|1918[ ). Given their 
ubiquity, astrophysical jets research has been a booming 
field of astronomy, especially after the advent of high- 
resolution imaging instruments and techniques (e.g., very 
long baseline interferometry, VLBI). Although observed in 
the minority of active galaxies (~ 15%; Krolik 1999), ex- 
tragalactic jets are some of the most pronounced morpho- 




Corresponding author: e-mail: mkarouzos@astro.snu.ac.kr 



Fig. 1 The active galactic nucleus 3C273 and its jet ob- 
served in three different wavelength regimes, the optical 
(left panel, HST), the X-ray (middle panel, Chandra), and 
the radio (right panel, MERLIN). Credit: NASA/STScI, 
NASA/CXC, MERLIN. 



logical features in AGN research. Their, presumably direct, 
connection to the active core and the SMBH residing there, 
makes them invaluable tools in the effort to characterize the 
properties and the underlying physics of activity in galax- 
ies. Although observable at different wavelengths (see Fig. 
[TJ, VLBI observations enabled the direct imaging of AGN 
jets at the highest possible resolutions and thus the detailed 
study of their properties. 

VLBI techniques have allowed the mapping of jets in 
the radio. A result of this was the discovery of the knotty na- 
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Fig. 2 An illustration of the AGN core, comprising the 
SMBH, the accretion disk, line emission regions, and the 
jet. Different jet regions and the effects observed in them 
are presented. Illustration reproduced from Lobanov (2007) 
(adapted from Marscher|2005| l. 



ture of AGN jets, observed as a consecutive series of bright- 
ness maxima and minima. Observations at different times 
revealed that these brightness maxima actually move. The 
investigation of extragalactic jet kinematics had begun. One 
of the most prominent discoveries, related to jet kinematics, 
was that of the superluminal movement of these jet compo- 
nents, which was first theoretically predicted by 
|Rees| ( fl966l l. The effect is a combination of relativistic ex- 
pansion speeds (close to the speed of light) and the projected 
geometry onto the plane of the sky. Superluminal movement 
was indeed detected, first i ndirectly ([Whitney et al.||1971 1, 
and then by direct imaging ( |Pearson & Readhead|1981 1. Su- 
perluminally moving components have become the staple of 
blazars, i.e., radio-loud AGN seen at the smallest viewing 
angles to their jet axis. 

Jet kinematics, as investigated through the study of dis- 
tinct components, is currently explained in terms of the 
shock-in-jet model (e.g., IMarscher & Gear||1985| also see 



Fig. |2j, where the observed jet knots are manifestations of 
shocks propagating at relativistic speeds down the jet. As 
discussed above, beaming and projection effects regulate the 
observed properties of the jets. Analysis of samples of ac- 
tive galaxies containing large number of sources, selected 
following stringent criteria, has been of fundamental impor- 
tance to t his end (e.g., |Ghiseilini||1993[ Vermeulen & Co- 
hen|1994||Vermeulen|1995[|Taylor et al.|1996||Hough et al. 
2002| ILister & Homan||2005[ |Britzen et al.||2007b| >. There 
has been a continuous effort to distinguish whether the dif- 
ferent types of radio-loud active galaxies (e.g., quasars, BL 
Lacs, etc.) and their jet properties are an effect of different 
viewing angles, or whether these objects have intrinsically 
different properties. The current paradigm is that the dif- 
ferences observed can be attributed to geometrical effects, 
although some indications to the contrary also exist (e.g., 
|Gabuzda|fT995j |Gabuzda et al.||2000l |Britzen et al.l|2009] 



Britzen et al. 2010b I. For example, using a sample of 39 su- 



perluminal sources Ghisellini et al. ( 1993 i find no apprecia 



ble difference between the Doppler factors distributions of 
BL Lacs and quasars, with radio galaxies showing smaller 
values. 



2 A different paradigm? 

Previous results concerning both individual sources as well 
as statistical samples have shown that flat-spectrum AGN 
predominantly show apparent superluminal motions, usu- 
ally interpreted in terms of radial, ballistic motions. Con- 
trary to this, recent detailed studies of a number of BL Lac 
objects suggest a set of somewhat different kinematic prop- 
erties (see Fig . [3]> . 

Through the (re-)analysis of more than 90 epochs of 
global VLBI and VLBA data of the BL Lac source S5 
1803+784, a new kinematic scheme for that source was re- 
vealed (Britzen et al. 2010a). All components in the inner 
part of the jet (up to 12 mas) appear to be quasi-stationary 
in terms of their core separation. This behavior is seen at all 
frequencies studied by the authors (1.6 - 15 GHz). In con- 
trast to this, the components show strong changes in their 
position angles, implying a dominant component of motion 
perpendicular to the jet axis. 

Britzen et al. (2010a) also studied the jet ridge line prop- 



erties of S5 1803+784. A jet ridge line at a given epoch is 
defined as the line that linearly connects all component po- 
sitions at that epoch (for some examples of jet ridge lines, 
see Fig. [5]). The jet ridge line of S5 1803+784 is found to 
change in an almost periodic manner, starting by resem- 
bling a straight line, evolving into a sinusoid-like pattern, 
and finally returning to its original linear pattern, although 
slightly displaced from its original position. The authors 
calculate a period of ~ 8.5 years for the evolution of the 
jet ridge line. 

Finally, the authors find that the jet changes its apparent 
width (in a range between a few and a few tens of degrees) 
in an almost periodic way, with a timescale similar to the 
one found from the evolution of the jet ridge line. All of 
the above properties support a new kinematic scheme for 
1803+784, where components follow oscillatory -like tra- 
jectories, with their movement predominantly happening 
perpendicular to the jet axis rather than along it. In addi- 
tion, the jet appears at times to form a wide channel of flow, 
while changing its width considerably across time. 

However, S5 1803+784 is not the only exception. The 
sources S5 0716+714 ( |Britzen et aT] [2009] l and PKS 



0735+178 ( |Britzen et al.|2010b| l, both classified as BL Lac 
objects, show similar behavior, with strong non-radial mo- 
tions discovered in their jets. Although not as pronounced as 
in S5 1803+784, but also lacking the unprecedented num- 
ber of VLBI epochs of S5 1803+784, both sources show jet 
components in their inner-jets that appear quasi-stationary 
with respect to the core. In addition, PKS 0735+178 appears 
to interchange between the classical blazar-like behavior of 
superluminally moving components and the behavior seen 
in S5 1803+784. PKS 0735+178 might be the key to our un- 
derstanding of these non-radial motion dominated jets and 
potentially implies a geometrical effect as the origin of the 
observed properties. 
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Fig. 3 A comparison of the kinematic properties of two fiat-spectrum radio-loud AGN, 1226+023 (left) and S5 1803+784 
(right). The plots show the core separation of the sources' jet components as a function of time. 1226+023 shows clear 
outward motion at superluminal speeds, while S5 1803+784 shows most components remaining at almost constant core- 
separations. 1226+023 plot reproduced from Kellermann et al.|2004 and S5 1803+784 from Britzen et al. 2010a For both 
plots, VLBA data at 15 GHz are used. 



2.1 Motivation 



Under the unification scheme of active galaxies (e.g., |An- 
1993] |Urry & Padovanil[T995) , both BL Lac ob- 



tonucci 



jects and fiat-spectrum radio-quasars (FSRQs) are believed 
to be active galaxies for which the viewing angle to their jet 
axis is very small, leading to strong relativistic effects. In 
contrast to the latter, BL Lac objects are additionally char- 
acterized by the absence of strong emission lines. Assum- 
ing similar viewing angle distributions, in the last several 
years the common classification of blazar has often been 
used to describe members of either class, also in terms of 
their jet properties. However, this phenomenological unifi- 
cation of FSRQs and BL Lacs comes into question in light 
of the recent investigation of sources like S5 1803+784, 
PKS 0735+178 and S5 0716+714. The question arises 
whether the peculiar kinematic behavior seen in these ob- 
jects is revealed due to the unprecedented richness of the 
datasets available, and is therefore relevant for all flat- 
spectrum radio- AGN, or whether BL Lacs are characterized 
by a genuinely different set of kinematic properties. 

We investigate the relevance of the above described, 
new, kinematic scheme for flat-spectrum radio- AGN, while 
investigating the apparent divide between BL Lac and 
FSRQ jet kinematics. We use the Caltech Jodrell-bank flat- 
spectrum (CJF) sample to statistically investigate and as- 
sess the similarity, or divergence, of the kinematic and mor- 
phological properties between the two distinct sub-samples 
of FSRQs and BL Lac objects in the CJF. We want to test 
whether jet components of BL Lac objects do show slower 
apparent speeds with respect to their cores compared to FS- 
RQs. Furthermore we are interested in the phenomenon ex- 
hibited in S5 1803+784 of a, at times, very wide jet, with a 
jet ridge line intermittently resembling a sinusoid-like mor- 



phology. To do that we use tools that extract information 
from the jet ridge line of the sources and characterize there- 
fore the jet globally. This allows for a investigation that is 
mostly independent of component modeling and component 
cross-identification. 

The paper is organized as follows: in Sect. [3] we de- 
scribe the CJF sample, in Sect. |4] we shortly describe the 
data and tools used, in Sect.[5]we present the analysis of our 
data and the results, and in Sect. [6] we discuss our results 
and give some conclusions. Throughout the paper, we as- 
sume the cosmological parameters Hq = 71 kms^'Mpc -1 , 
Q M = 0.27, and Q A = 0.73 (from the first-year WMAP 
observations; Sperg el et al.|2003| l. 



3 The CJF sample 



The CJF sample ( |Taylor et al.|1996) is a flux-limited radio- 
selected sample containing 293 radio-loud active galaxies 
selected from three different samples (see Table [T] and 
Britzen et al. 2007b ). The sources span a large redshift range 
(see Fig.|4|, the farthest object being at a redshift z = 3.889 
(1745 +624; [Hook et al.|19"95] l and the closest at z = 0.0108 
(1146+596; de Vaucoule urs|1991 1. The average redshift of 
the sample is z avg = 1.254, z B LLac, a v g = 0.546, z R c,av g = 
0.554, and ZFSRQ.avg - 1-489 for BL Lacs, radio galax- 
ies, and FSRQs, respectively. All the objects have been ob- 
served with the VLBA and/or the global VLBI network. 
Each source has at least 3 epochs of observations (with a 
maximum of 5 epochs) and has been imaged and studied 



kinematically (|Britzen et al.||1999| |Britzen et IT] [2007b 



Britzen et al.| |2008). The X-ray and gamma-ray properties 



have been studied and correlated with their VLBI proper- 
ties ( Britzen et al.|2007a and Karouzos et al.|20TT respec- 
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Fig. 4 Redshift distribution of radio galaxies (grey 
blocks), BL Lacs (dashed line), and FSRQs (solid line) in 
the CJF sample. 



tively). The evolution of active galaxies, in the context of the 
merger-driven evolution scheme (e.g., Hop kins et al.|2 006 ), 
has also been investigated with the help of the CJF, identify- 
ing candidate CJF sources in different evolutionary stages, 
including new binary black hole candidates (Karouzos et al. 
|20T0l >. 



Table 1 CJF sample and its properties. 



Frequency (MHz) 


4850 


Flux lower limit @5 GHz 


350mJy 


Spectral Index 


0,4850 > _q 
"1400 - " 


Declination 


S> 35° 


Galactic latitude 


\b\ > 10° 


# Quasars 


198 


# BL Lac 


32 


# Radio Galaxies 


52 


# Unclassified 


11 


# Total 


293 



4 Data 

The work presented here is heavily based on the kinematic 



analysis of the CJF sample (Britzen et al. 2007b Britzen 
|et al.|2008] >. Five CJF sources were initially excluded from 
any further analysis due to problematic observations (0256+ 
424, 0344+405, 0424+670, 0945+664, 1545+497) and 
therefore are not considered here. In total, 288 sources are 
considered and analyzed in the following sections. Of these, 
according to the optical classification from |Britzen et al.| 
P007b| l, 196 are classified as quasars, 49 as radio galaxies, 
33 as BL Lac objects , and 10 are not classified. 

Due to the scope of the CJF program, the identification 
and analysis of pc-scale jet component kinematics has fo- 
cused on the part of the jet that is beamed towards us. For 
a number of sources, several components belonging to the 
counter-jet have been identified. However, for these sources 
cross-identification of the counter-jet components over 



epochs has not been carried out. For this reason, and given 
the nature of the analysis that we undertook (see below), we 
have excluded all counter-jet components in the following 
investigation. 

The tools that we use for the analysis of the CJF jet ridge 
lines, in the context described in Sect. [2] are described be- 
low: 

- Monotonicity Index, M.I. 

- Apparent Jet Width, dP 

- Apparent Jet Linear Evolution, A£ 

We note that all the above measures, as their names imply, 
refer to values projected onto the plane of the sky. For the 
following sections, we adopt the basic assumption that BL 
Lacs and FSRQs are seen at the smallest viewing angles. As 
we are interested in the comparison between FSRQs and BL 
Lacs, deprojection of the jet properties investigated here is 
not critical. For the following analysis, and for the sake of 
brevity, we shall drop the characterization of "apparent" for 
each of these values, although this will be implied through- 
out. 



4.1 Monotonicity Index, M.I. 

Detailed VLBI investigations in the past have shown that a 
large number of active galaxies exhibit bent or otherwise 
non-linear jet morphologies on various scales. Individual 
sources like S5 1803+784 and PKS 0735+178, as well as 
others (e.g. 3C 345;|Lobanov & Zensusp999l B0605-085; 
Kudryavtseva et al.|2010| i, have been extensively studied to 
understand the origin of this bending. In this context, we 



are interested in quantifying the bending of the whole jet 
ridge line. Moreover, we want to differentiate between a 
"monotonically-bent" jet, i.e., a jet that is bent in only one 
direction (see Fig. [5] right), as opposed to a more sinusoid- 
like morphology (similar to what is seen for some epochs of 
S5 1803+784, see Fig.[5} left). This is done through the use 
of the monotonicity index, M.I.. 

We quantify a sinusoid-like morphology of a jet by iden- 
tifying the local extrema in a given jet ridge line, in the core 
separation - position angle plane. For an epoch i, a compo- 
nent m exhibits a local extremum under the definition 

6m;extr '■ \@m ~ &m±l\ S I0(dd m + d6 m ±\), (1) 

where 6 m and d9 m denote the position angle of component 
m and its uncertainty. 9 m is calculated as 
X 



6 m = arctan 



Y' 



(2) 



where X and Y are cartesian coordinates on the plane of the 
skjQ Having calculated the number of extrema for a given 
jet ridge line at an epoch i, we define the M.I. as 

number of extrema 



M.I. = 



N -I 



(3) 



1 We assume a conservative value of 10tr in Eq. ^ (where a here is 
defined as the sum of the position angle errors for 9 m and 8 m± \ ) as the lower 
limit for flagging an extremum. This reflects an original underestimation 
of the individual 6„, uncertainties in the component fitting of Britz en et al.| 
j2007b| . 
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Fig. 5 Radio maps at 5 GHz of 1928+738 (left) and 0700+470 (right), with jet ridge lines superimposed. Two examples 
of a sinusoid-like jet morphology (with M.I. =0.75, left) and a single-bent one (with M.I.=0, right). VLBI maps from 
IBritzen et alT| < |2007b) . 



where N is the total number of components at that given 
epoch. This is a crude calculation, but it can give us a mea- 
sure of how the bending of the jet behaves along the jet. For 
M.I. values close to one, the jet more closely resembles a 
sinusoid. An M.I. value close to zero reveals a monotonic, 
single-bend, jet morphology. We normalize for the number 
of components N to account for longer, or shorter, jets and 
to enable comparison between different sources. An exam- 
ple of such differences in the jet bending is shown in Fig. 

s 

4.2 Apparent jet width, dP 

For an epoch i and a jet consisting of N components char- 
acterized by their distance from the core and position angle 
(r,, Oi), we identify the components with the maximum and 
the minimum position angles. The apparent width of the jet 
dP, measured in degrees, is then calculated as (see Fig.|6]i: 

dPt = ff; ,ax - 0™\ (4) 

while the error is calculated by the propagation of errors 
formula. One obvious drawback for the above definition is 
the non-localized nature of this measure. Using two differ- 
ent components, at different core separations, gives us only 
an approximate notion of the width of the flow. By using the 
definition above, we seek to quantify the opening of the jet 
flow and identify the effect of a channel-like jet, as seen in 
the case of S5 1803+784. 

4.3 Apparent jet linear evolution, At 

The linear evolution across all available epochs and for all 
components of a jet can be calculated, ultimately produc- 
ing a value that reflects the total linear displacement of the 
whole jet ridge line. We use plane-of-the-sky coordinates 
(Xj, Yj) to calculate the linear displacement of component m 
between epochs i and 

If = VAX 2 + AY 2 . (5) 



Clean map. Array: BFHKLMNOPS . 
0700+470 at 4.991 GHz 1995 Sep 04 ClP = 0max "©min 




Right Ascension (mas} 

Fig. 6 Radio map of 0700+470 at 5 GHz, exhibiting the 
calculation of the jet ridge line width. The two arrows de- 
note the components at maximum and minimum position 
angle, while the area in between is defined as the opening, 
or apparent width, of the jet ridge line. The VLBI core is 
found at (0,0) coordinates. VLBI map from |Britzen et al. 
( |2007bD 

To calculate the total displacement of the whole jet ridge 
line we then need to sum up over all components and all 
available observing epochs: 

i-l m 

We need to account for both the different time span of ob- 
servations, as well as the different number of components. 
We therefore define the jet linear evolution Al, measured in 
parsecs per unit time and per component, as follows: 

Aii— 1 2m I 

At = = , (7) 

N(Ti-Ti) NdT 

where N here is the total number of components used across 
all the epochs. For the calculation of the displacement of an 
individual component I between two consecutive epochs, 
and consequently of Al, the cross-identification of compo- 
nents across epochs is necessary. Therefore, for Al, unlike 
the previously discussed jet ridge line characteristic values, 
the actual identification of components is important. 
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A( essentially reflects the apparent speed distribution of 
all cross-identified components of the jet and therefore rep- 
resents a value characteristic for the whole jet, rather than 
for any individual component. By summing up all compo- 
nents and epochs we trade temporal and positional resolu- 
tion for a universal treatment of the entire jet. In this way 
we can test whether the kinematics of BL Lac objects is 
fundamentally different than that of FSRQs while averag- 
ing out localized properties of individual components. A£ 
can be seen as a mean jet component speed, with the differ- 
ence that it is acquired through averaging not only over all 
components, but additionally over all available epochs. Al- 
though it certainly reflects a measure of the outward motion 
in BL Lac jets, the calculation of Al is done in such a way 
that the potential curvature of the components' trajectories 
is taken into account. This separates Al from a simple linear 
regression fit to the core-separation versus time diagrams 
usually employed to calculate outward velocities, making it 
sensitive to non-radial motions, that are otherwise missed. 



5 Analysis and results 

Since the CJF sample spans a large redshift range, com- 
bined with the fact that we are using apparent quantities, 
i.e., measures that are projected on the sky, different linear 
scales will be probed at different redshift ranges. To avoid 
this effect we impose two constraints to our data. 

- As most CJF sources have available spectroscopic red- 
shifts, we can calculate their projected jet length in lin- 
ear scales. We can then impose a constraint to the max- 
imum core separation out to which a component is con- 
sidered in this investigation. Thus we limit ourselves to 
the same linear scales for all of the CJF objects, albeit 
with different linear resolution. As we are interested in 
the behavior of BL Lacs, we use them as the basis of 
our decision. For S5 1803+784, peculiar kinematics of 
the components is observed for the inner-most part of 
its jet, out to 6 mas. For this source's redshift (z=0.68), 
this translates roughly to 41 pc. Similar to S5 1803+784, 
PKS 0735+178 and S5 0716+714 also show the most 
prominent evolution of their jet ridge line out to ~ 6 mas 
(39.6 pc and 33 pc, for their respective redshifts). Given 
that most of the BL Lacs in our sample are at z < 1, 
the length scale implied by the three sources mentioned 
above (the only objects whose jet ridge line has been in- 
vestigated in detail) ensures that the inner-jet of sources 
at z < 1 is covered. We adopt therefore an apparent core 
separation limit of 40 parsecs. 

- As we are interested in the behavior of BL Lac objects, 
we can limit our comparisons to the redshift range [0,1]. 
In this way we can partially account for the different lin- 
ear scales, as well as avoiding any possible evolution of 
the values we are investigating as a function of redshift. 
There are 21 BL Lac sources with available redshift in- 



formation, with 17 of them being at z < 1. There are 56 
FSRQs in the same redshift range. 

5.1 Jet ridge line morphology 

We investigate how, if at all, a jet resembles a sinusoid (as in 
the case of S5 1803+784). In order to do this we calculate 
the M.I. from the jet ridge lines of our sources. The M.I. is 
calculated for the whole jet (we do not apply the core sepa- 
ration limit of 40 pc), as we are interested in the morphology 
of the entire jet ridge line, rather than a localized property. 
Additionally, a large number of components is needed for a 
robust interpretation of the M.I., we therefore take into ac- 
count only epochs with at least three components identified. 
Finally, we define the maximum M.I. value for each source 
across the available epochs. Fluctuations of the M.I. of a 
source between epochs are low. 

Under the above constraints, we find that there are in to- 
tal 46 sources with M.I. > 0.5. This translates to 37% of the 
total number of sources used here. BL Lacs show more often 
jets with sinusoid-like morphologies, compared to FSRQs. 
More than half of the FSRQs (~ 51%) have an M.I. value 
of 0, while two thirds (~ 66%) have an M.I. value lower 
than 0.5. In contrast to this, only one third of the BL Lacs 
(~ 33%) have M.I.=0 and less than half (~ 44%) of them 
have M.I.< 0.5. Figure|7]shows the histogram of the M.I. for 
the different types of objects. BL Lacs show their maximum 
in the [0.4,0.6) bin. In contrast, FSRQs have their maximum 
in the [0,0.2) bin. When considering the median values for 
the two classes of objects, BL Lacs and FSRQs differ at 
> 3cr level (0.500+0.014 compared to 0.000+0.004, respec- 
tively). It should be noted that a robust quantitative treat- 
ment of the M.I. is problematic, given the severe restrictions 
imposed by the small number of components per jet, as well 
as the limited temporal resolution of the observations used 
here. 
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Fig. 7 Monotonicity Index (M.I.) distribution for BL Lacs 
(filled line) and FSRQs (open line). The histogram has been 
normalized to area unity. 
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BL Lacs appear to more often show sinusoid-like curved 
jets, as seen for S5 1803+784 and 0716+714. 



that 47.3% of BL Lacs have jets wider than 20 degrees, as 
opposed to only 13.6% for FSRQs. 



5.2 Apparent jet width, P, statistics 

We calculate non-zero jet widths for a total of 559 epochs 
(multiple epochs per source). For one epoch of the source S5 
1803+784 and one epoch of the source 2200+420, a width 
of ~ 180 degrees is found. This is probably due to a problem 
in the original component fitting and we therefore exclude 
these two epochs in this analysis. Of the remaining 557, 277 
are larger than 10 degrees, while 125 are larger than 20 de- 
grees. This translates to 22.4% of the sample possessing jets 
with widths exceeding 20 degrees. 

We calculate an average jet ridge line width for our 
sources. We find that BL Lacs show significantly wider jet 
ridge lines (average value of 19.3° ± 1.8°), compared to FS- 
RQs (12.0° ± 0.4o; 4cr difference). A Student's t-test gives a 
significance >99.99% that BL Lacs and FSRQs show dif- 
ferent mean values. As discussed above, we also investi- 
gate sources in the redshift bin [0, 1 ] . In Table [2] we show 
the characteristic statistical parameters for the CJF jet ridge 
line widths for all sources with available redshift, as well as 
sources in the redshift bin [0,1]. 



Table 2 Characteristic statistical values concerning the 
width of the jet ridge line of sources with measured red- 
shifts. We give average and median values with uncertain- 
ties for FSRQs and BL Lacs. 



Jet Width C) 




All 




< z 


< 1 


Types 


FSRQ 


BL 


FSRQ 


BL 


# 


343 


63 


91 


58 


Average 


12.0 


19.3 


11.1 


20.4 


Error 


0.4 


1.8 


0.8 


2.0 


Median 


8.7 


12.5 


7.5 


13.8 


Error 


0.5 


2.2 


1.0 


2.5 



Focusing on the redshift bin [0,1], from Table [2] we see 
that BL Lac jet ridge lines appear significantly wider than 
their FSRQ counterparts (20.4° + 2.0°, compared to 1 1.1° + 
0.8°; > 4cr). The same behavior is seen when looking at the 
median values (albeit at the 2<x level). 

In Fig.[8]we show the distribution of jet ridge line widths 
for BL Lacs and FSRQs, in the redshift bin [0,1]. Both 
classes show similar distributions, with their maxima sit- 
uated at around 10 degrees. The distribution of BL Lac ap- 
pears to be wider and extending to larger widths, compared 
to FSRQs. FSRQs show a rather more contained distribu- 
tion to lower width values than BL Lacs. For the 125 epochs 
where CJF sources show apparent jet widths > 20 degrees, 
we find 14 BL Lacs and 27 FSRQs. When taking into ac- 
count the total number of each type of object, we calculate 
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Fig. 8 Jet ridge line width distribution for BL Lacs (filled 
line) and FSRQs (open line). For this histogram sources 
with redshifts in the [0,1] bin are used. 



We apply the two-sample Kolmogorov-Smirnoff (K-S) 
test to our data to see whether the BL Lacs and FSRQs are 
indeed different in their jet ridge line widths. Comparing the 
jet ridge line widths of BL Lacs and FSRQs (no redshift bin- 
ning) the K-S test gives a 0.7 x 10~ 3 % probability that these 
two sub-samples originate from the same parent distribu- 
tion. We also compare the sub-samples of BL Lacs and FS- 
RQs in the redshift bin [0,1]. In this case the K-S test gives a 
probability of 1% that BL Lacs and FSRQs are drawn from 
the same parent sample. We conclude that BL Lacs show 
a significantly different distribution of jet ridge line widths 
compared to FSRQs. In both average and median values, 
BL Lacs show substantially wider jets than FSRQs. 

5.3 Apparent jet linear evolution, A£, statistics 

Focusing on the final measure of the jet ridge line kinemat- 
ics, as studied in this paper, we want to investigate whether 
the stationarity of components, as observed in the case of S5 
1803+784, is commonplace among other BL Lacs. We use 
the linear evolution measure, as described in Sect. |4] to do 
this. We follow the same procedure as previously to check 
the statistics of the individual classes. 

In Table [3] we give the statistical properties of the total 
linear evolution of the jet ridge line (measured in parsecs 
per unit time and component) distributions for FSRQs and 
BL Lacs. On average, BL Lacs show weaker evolution of 
their jet ridge lines compared to FSRQs (~ 2<x difference). 
Looking at the median values, BL Lacs show the least evo- 
lution compared to FSRQs (0.24 + 0.06 and 0.414 + 0.020 
pc/yr/comp respectively; ~ 3<x difference). We also study 
the statistics of the sub-sample of CJF sources in the [0,1] 
redshift bin (also in Table[3]l. The behavior remains the same 
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as before, with differences between BL Lacs and FSRQs be- 
coming more pronounced (with a ~ 4<x difference in median 
values). 



Table 3 Statistical properties of the jet linear evolution 
distributions for FSRQs and BL Lacs. Average and median 
values are calculated for both all sources, as well as for 
sources in the redshift bin [0,1]. Only the inner part of the 
jet is considered (< 40 pc). 



Jet linear evolution {pc/yr/comp) 
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z 


< 


z< 1 


Types 


FSRQ 


BL 


FSRQ 


BL 


# 


171 


25 


44 


21 


Average 


0.470 


0.38 


0.51 


0.38 


Error 


0.016 


0.05 


0.04 


0.05 


Median 


0.414 


0.24 


0.50 


0.24 


Error 


0.020 


0.06 


0.04 


0.06 




distance to core 



Fig. 10 A sketch of the new paradigm for jet kinematics 
in BL Lac objects, compared to the established scheme of 
superluminal motion in blazars. BL Lac jet components ap- 
pear to move less outwards, while exhibiting stronger non- 
radial motion components. 



We now turn to the actual distribution of the total jet lin- 
ear evolution for the three types of objects, shown in Fig. [9] 
We once again only take into account sources in the redshift 
bin [0,1]. Quasars show a pronounced maximum around 
0.63 pc/yr/comp, extending out to 2 pc/yr/comp. In con- 
trast, BL Lacs show their distribution maxima around 0.25 
pc/yr/comp. Quasars appear to have a wider distribution of 
total linear evolution values, also showing the highest maxi- 
mum values among all three types (see Table[3]l. Once again 
we employ the K-S test to compare the distributions. For the 
whole sample (independent of redshift constraints), the K- 

5 test gives a probability of 4.3% that BL Lacs and FSRQs 
are drawn from the same parent population. Focusing on the 
< z < 1 sub-sample, we get a marginal 5.4% probability 
that BL Lacs and FSRQs stem from the same parent sample. 
These results give a positive answer as to whether BL Lacs 
show less absolute linear evolution of their jets compared to 
their FSRQ counterparts. 

6 Discussion 



Figure 10 shows a sketch of the new paradigm as was origi- 
nally seen in the BL Lac object S5 1803+784 and the effects 
of which we tried to investigate in the previous sections. 
Our statistical investigation, part of which is presented here, 
lends support to a different kinematic scheme for BL Lac 
objects. Indeed we find that BL Lac jet components show 
significant non-radial components in their motion, 
compared to those of FSRQs. It should however be noted 
that there is a substantial part of the total CJF sources that 
show similar, "BL Lac -like" behavior. This implies the pres- 
ence of a more universal effect than simply a property char- 
acteristic to BL Lacs. 



6.1 Possible explanations 

It is of great interest to try and explain the kinematic be- 
havior seen in a large number of the CJF sources. Combin- 
ing the percentages calculated from both apparent width and 
apparent linear evolution we get an approximately 30% of 
the sample showing wide jets that show slow outwards mo- 
tions. There also seems to exist an apparent preference for 
BL Lacs to show this behavior over FSRQs. One obvious 
factor that should play a deciding role for the kinematics of 
jet components is the viewing angle under which a source is 
observed. BL Lacs and FSRQs are assumed to be seen jet- 
on, at the smallest viewing angles, with steeper-spectrum 
quasars and radio galaxies being sources observed at pro- 
gressively larger viewing angles (e.g., Antonucci 1993 



Urry & Padovani||1995|l. Could a viewing angle difference 



explain the kinematic differences seen in some of the CJF 
sources and in particular observed between BL Lacs and 
FSRQs? Assuming that jet components follow ballistic, lin- 
ear paths, one can expect that viewed at smaller viewing 
angles (smaller than the critical angle l/y, with y being the 
Lorentz factor of the flow), the components are observed to 
cover smaller distances than if seen edge-on. Therefore, the 
slower speeds, and hence smaller total linear evolution of 
their jet ridge lines could be explained in terms of a sys- 
tematically smaller viewing angle for BL Lacs compared to 
FSRQs. 



However, Hovatta et al. ( 2009 ) used long timescale vari- 
ability of AGN, together with jet kinematics, to derive 
Doppler factors and consequently viewing angles of a sam- 
ple of 87 AGN. They found that FSRQs actually show 
smaller mean viewing angles than BL Lacs (a result also 
found by Lahteenmaki & Valtaoja 1999). Although there 
are certain limitations to the method used by the authors, 
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Fig. 9 Distribution of the total linear evolution of the CJF jet ridge lines for FSRQs (open line) and BL Lacs (filled line). 
We use sources with redshift lower than 1 . 



their results are difficult to reconcile with a viewing-angle- 
dependent explanation of our results. 

Our results would therefore indicate the need for an ad- 
ditional, potentially geometric, effect in play. This would be 
in agreement with that low-speed components in FSRQ and 
BL Lac jets should appear so because their pattern Lorentz 
factor is lower than the bulk Lorentz factor of the jet (e.g., 
Cohen et al. 2007). Alternatively, it should be considered 
whether there is some systematic bias in the way Doppler 
factors are estimated that would lead to an over- or under- 
estimation of the viewing angles of BL Lacs and FSRQs, 
respectively. 

The large jet widths are more difficult to explain. An 
additional mechanism or effect must be introduced to pro- 
duce a wide jet. Precession of the jet axis, or assuming that 
the components follow non-ballistic and non-linear paths, 
can lead to such jet properties (e.g., Steffen et al.|[T995b ; 



Gong|2008| |Roland et al.|2008l |Gong et al.|20TT) 7ln these 
models, jet components follow highly curved trajectories, 
which, viewed face-on, would give the impression of a 
wider distribution of jet component position angles. A pre- 
cessing jet, or rather a precessing jet nozzle, would imply 
that different components follow different trajectories. That 
would result in a changing jet component position angle dis- 
tribution, more pronounced at smaller viewing angles. 

That more BL Lacs show M.I. closer to one compared to 
FSRQs, implying a highly curved, sinusoid-like ridge-line, 
lends additional support to this scenario. That leads us to the 
conclusion that either (1) there is a sub-set of FSRQs and BL 
Lacs (~ 30% in our sample) for which non-ballistic and/or 
precession effects play an important role, or more gener- 
ally that (2) for these sources the viewing angles are below 
some critical angle 6 C that allow us to witness and thus char- 
acterize the helical or "non-ballistic" structure common in 
all radio- AGN jets, but which would otherwise be blended 
out by the projection effects at viewing angles > C . 



A final scenario that needs to be discussed is whether 
a systematic difference in the Lorentz factors y between 
FSRQs and BL Lacs could produce the effects observed 
here. It has been argued that BL Lacs show smaller Lorentz 
factors than FSRQs (e.g., |Morganti et al.| |1995 



Urry & PadovamfT995| |Hovatta et al.||2009| >. If this is true 



(and not a selection bias in the samples used) then that 
would naturally explain the slower components in BL Lacs. 
The above in turn ties in to the currently accepted unifica- 



tion scheme (e.g., Urry & Padovani 1995 1. In that paradigm, 
BL Lacs and FSRQs are drawn from two, presumably differ- 
ent, parent samples of Fanaroff-Riley I (low luminosity) and 
II (high luminosity) galaxies (FR; FanarofF & Riley||1974 



[ Padovani &TJrryl[T99T] [Capetti & Celotti||1999| |Xu et al.| 
2009), respectively. If that is true, the differences seen in the 
jet ridge line properties of BL Lacs and FSRQs should then 
translate to differences between FRI and FRII jet kinemat- 
ics. Interestingly, studies of FRI and FRII kinematics have 
shown that both types have similar parcec-scale jet Lorentz 
factors (e.g., Giovannini et al.|2001 1, contradictory to what 
we find here, as well as in other studies mentioned above, 
concerning the parsec-scale jet speeds in BL Lacs and FS- 
RQs. 

6.1.1 Helical jets 

There is a growing number of sources for which helical 
models are employed to explain their kinematic, flux, and 
evolution properties (e.g., 3C 273, |Lobanov & Zensus 
fRoland et al.||2008 



1803+784, 



3C454.3, pan et al. 



2001 



2007 



0605-085, |Kudryavtseva etaL]|2010| 1632+382, pu et al. 
2010 etc.). A helical, bent, trajectory has been employed 
for 3C345, showing that a constant Lorentz factor of about 
10 can explain the observations. Apparent acceleration of 
the components is an effect of the change of the viewing an- 
gle ( |Zensus et al.|1995| ). Further evidence supporting an in- 
trinsic helical structure for AGN jets (as traced by the mag- 
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netic fields threading them) is given by polarization stud- 
ies, and in particular Faraday rotation gradients and circular 
polarization found in several sources (e.g., Gabuzda 2005; 
IGabuzda et al.|2008) . 

We can use a simple helical jet model and try to see 
whether a helical geometry, combined with projection ef- 
fects, can account for the results of the statistical analysis 
above, namely that BL Lacs are expected to have wider jets, 
than FQSOs, and also show weaker total linear evolution of 
their jets. 

The helical model of Steffen et al. (1995b), applied in 
that case to the QSO 3C345 (the same model has also been 



applied to the BL Lac source S5 1803+784, Steffen et aL| 
[199311 QSO B3 1633+382, |Liu et al.||201oj can be en> 
ployed to that end. The motion of a component down the jet 
is determined by the conservation of a number of quantities 
(e.g., specific kinetic energy, specific angular momentum, 
etc.). In particular, [Steffen et al7| (|T995b > find that the best fit 
to 3C345 provides a case where the specific kinetic energy 
of the component, the specific angular momentum, as well 
as the opening half-angle of the jet are conserved (Case 3 in 
that paper). The geometry of the model can be seen in Fig. 

El 



VLBI-core 
position 




helical path 



Fig. 11 The geometry of the helical jet model of Steffen 



et al. ( 1995b i. The VLBI jet core is put at the origin of the 
cylindrical coordinates system (r, 0, z). The jet axis (parallel 
to the z-axis) forms an angle 9 to the line of sight. Figure 
reproduced from Steffen et al. ( 1995b) >. 



In this context, the trajectory of a single component in 
the jet can be defined by the following set of parameters: the 
Lorentz factor and angular momentum of the component, 
the opening half-angle of the jet iff, and the inclination angle 
of the jet axis 9. Of particular interest for this investigation 
is the opening half-angle of the jet tfr, which defines how 
helical or linear the trajectory of a component is and the 
inclination angle of the jet axis, or viewing angle, 9, which 



defines the projection effects at play. In Fig. 12 examples of 



component trajectories (in sky coordinates) are given for 6 
different values of i//. For each tfr we allow for three different 
values of the viewing angle 9. Both helicity and projection 
effects can be seen. A relativistic component with a speed 
u = 0.995c is assumed. 



We can use this model to study whether and how the 
projected width of a jet depends on its helicity and the angle 
at which it is viewed. By definition, this model calculates 
the evolution of a single component. This is of course not 
the case of what is observed for most sources. The width 
of the CJF jet ridge lines is calculated as the difference (in 
degrees) between the components with the maximum and 
minimum position angle at a given epoch. Making the sim- 
plifying - and potentially simplistic - assumption that all 
components follow the same path (therefore implicitly as- 
suming that i/r, 9, and u do not change between different 
components), this is equivalent to finding the maximum and 
minimum position angle for one component across different 
epochs. 

In Fig. [13] (left) the evolution of the position angle of a 
component is shown as a function of time, for a helical jet 
(ijj - \) and for 4 different viewing angles. It is obvious that 
for smaller viewing angles, the modulations of the position 
angle become stronger, resulting in substantially wider po- 
sition angle distributions than for the case of large viewing 
angles. 

We can calculate the projected width of each loop of 
the helix. The maximum value of these widths (for a given 
viewing angle) is representative of that helix and therefore 
it is interesting to investigate its behavior as a function of 



9 and iff. This is shown in Fig. 13 (right). As expected, for 
highly helical jets (small values of i^), smaller viewing an- 
gles produce substantially wider jets. The highest value of 
apparent width at the smallest viewing angle is produced by 
the if/ — 1 jet, while for slightly larger viewing angles (up to 
~ 15 degrees) the iff = 2 jet appears wider. For tfr > 2 zero 
apparent jet width values are not calculated. This is because 
of the method used to calculate the apparent width and re- 
flects the fact that for these jets no single loop is produced in 
the rest-frame of the jet (i.e., these jets do not have a helical 
structure). 

In conclusion, helical jets appear to reproduce well the 
effects seen for the width of the CJF jet ridge lines. A com- 
bination of a highly bent trajectory of the components, to- 
gether with the effects of the projected geometry, can ex- 
plain the extremely wide jets seen for BL Lacs (as it can be 



seen in Fig. 13 right). In this context, the subset of sources 
for which significant apparent jet widths are seen (i.e., dP 
> 20 deg) are the sources at the smallest viewing angles. 
The trend for BL Lacs to show apparently wider jets than 
FSRQs would then imply that BL Lacs should be at system- 
atically smaller viewing angles than the FSRQ sample. 

The origin of such a helically structured jet is an all to- 
gether different question. Such helical structures/trajectories 
have been explained by a number of different models. An of- 
ten employed concept is that of a precessing jet, combined 
with the non-ballistic motion of ejected components. For a 
detailed discussion see |Britzen et al.| ( f2010a[ ). 
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Fig. 12 Examples of different component trajectories for different opening half-angles: (a) = 0.2, (b) 1, (c) 2, (d) 5, 
(e) 9, and (f) 14 degrees, as well as three different inclination angles 8 of the jet axis: (red) 2, (black) 45, and (green) 85 
degrees. X and Y are sky coordinates, the plane perpendicular to the observer's line of sight. 




Fig. 13 Left: Position angle as a function of time for a helical jet (if/ = 1) and for 4 different viewing angles: (red) 2, 
(blue) 10, (orange) 45, and (green) 85 degrees. Right: Apparent jet width as a function of the viewing angle, for different 
opening half-angles: (red) 0.2, (blue) 1, (green) 2, and (orange) 5. 



6.2 Comparison with previous studies 

It should be noted that this is the first time the jet ridge lines 
of a sample of AGN this size have been explicitly and ex- 
clusively studied in a statistical manner. While the ridge line 
of a jet is not a new concept, so far only the jet ridge lines 
of individual sources (both galactic and stellar) have been 
studied (e.g.,|Condon & Mitchell| 19841 |Steffen et al.|1995b 



ILister et al.|20031|Lobanov et al.|20'06| |Britzen et al.|2009 



|Perucho et al 120091 pstereTa l. 2009 ; Bntzeii et al.|2010a| ). 



The definition of the jet ridge line can differ from study to 
study. We define the ridge line of a jet at a certain epoch 
as the line that linearly connects the projected positions of 
all components at that epoch. Alternatively, one can use an 
algorithm to more directly extract the jet ridge line from the 
VLBI map of each epoch (e.g., Pushka rev et al.|2009 1. 



The one obvious direct comparison that can be drawn is 



to the MOJAVE / 2cm sample (Lister & Homan 2005 ; Paper 
I). As the work presented here is prototype in its concep- 
tion, it is difficult to draw direct comparisons to the MO- 
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JAVE sample. The CJF sources appear to be significantly 
slower than what is seen in the MOJAVE sample (e.g., com- 
pare Fig. 7 from |Britzen et al . 2008 and Fig. 7 from |Lis^| 
ter et al. 2009). In the distribution of apparent speeds, for 



the CJF there is a maximum at around 4c and then a turn 
down, with very few sources above 15c. In contrast, the dis- 
tribution for the MOJAVE sources shows a maximum at 10c 
with a fair number of sources showing speeds greater than 
20c. This can be explained in terms of the higher tempo- 
ral resolution of the MOJAVE sample (owning to the larger 
number of epochs per source) that might allow the detec- 
tion of such faster motions. It is interesting to note that no 
clear distinction has been made in the MOJAVE sample be- 
tween BL Lacs and FSRQs, probably due to the relatively 
small number of BL Lacs. [Lister e t al. (2009) however do 
mention that BL Lacs are more often found to exhibit "low- 
pattern speeds", i.e., components with slow motions, sig- 
nificantly smaller than others in the same jet. Similar evi- 
dence has been previously found in different VLBI samples 



(e.g., Jorstad et al.|2001 ; Kellermann et al.|2004|). Combined 



with the argument that BL Lac jets exhibit lower Lorentz 
factors, as a result of a possible correlation between intrin- 
sic AGN luminosity and jet Lorentz factor (e.g., Morganti 



|et al.|1995] ), the above are in agreement with CJF BL Lacs 
showing on average weaker linear evolution than FSRQs. 
It should be noted that the agreement with previous studies 
concerning the outward motions of BL Lacs objects sup- 
ports the robustness and reliability of our method. 

A final note concerning the comparison between our re- 
sults and these of the MOJAVE sample pertains to the ob- 
servations frequency of each sample. Unlike the CJF, the 
MOJAVE uses VLBI observations at 15 GHz, meaning that 
a potentially different regime of the jet is probed, compared 
to what is studied in the CJF. Two-zone models, such as the 
"two-fluid" model (e.g., |Sol et al.|1989||Pelletier & Roland 
1989), predict a fast jet spine, embedded in a slower mov- 
ing sheath. This would imply that the MOJAVE potentially 
probes deeper, thus faster, jet layers compared to what is 
seen in the CJF sample. Moreover, higher observation fre- 
quency also translates to higher angular resolution and 
therefore to probing effects closer to the core and also en- 
abling the detection of faster moving components. In this 
case, it would be interesting to apply the methods used in 
this paper for the MOJAVE dataset, as it is expected that the 
curvature of the jet (and therefore potentially the angular 
properties and evolution of the jet ridge line) should play an 
increasingly important role at closer core separations. 

7 Conclusions 

We develop a number of tools to investigate both the mor- 
phology and kinematics of the CJF jet ridge lines: 

1 . Monotonicity Index 

2. Apparent jet width 

3. Apparent jet linear evolution 
We summarize our results: 



- Using the M.I., we find that BL Lacs jet ridge lines more 
often resemble a sinusoidal curve compared to FSRQs. 
In contrast, 2/3 of the FSRQs have an M.I. lower than 
0.5 (indicating fairly monotonic jets). 

- 22.4% of the CJF sources have apparent jet widths larger 
than 20 degrees. 47.3% of BL Lacs, over 13.6% for FS- 
RQs, have dP > 20 degrees. 

- BL Lacs exhibit substantially apparently wider jets than 
FSRQs. This effect persists for a smaller redshift range. 
This supports the effects seen in individual BL Lac ob- 
jects (i.e., 1803+784, 0716+714, etc.). 

- The distribution of apparent jet ridge line widths for 
BL Lacs appears to extend towards higher values, with 
FSRQs mainly contained at lower values. A K-S test in- 
dicates a 0.7 x 10 3 % probability that FSRQs and BL 
Lacs are drawn from a single parent population. 

- BL Lac objects, on average, show weaker apparent lin- 
ear evolution of their jet ridge lines compared to FSRQs. 

- We use a helical jet model to show that a helical geom- 
etry combined with a small viewing angle can explain 
the large widths seen in BL Lac objects. 

In conclusion, by statistically analyzing the CJF sam- 
ple, we provide detailed insight concerning the morphology 
and evolution of AGN jet ridge lines. The statistical investi- 
gation of the CJF sources lends independent support to the 
different kinematic scenario recently seen in a number of 
BL Lac objects (180 3+784, |Britzen et al.|2010a| 0735+178, 
Britzen etal.|2010b| ). 25%-30% of the CJF sample show ap- 
parently considerably wide jets and strong apparent width 
evolution. BL Lac objects appear to deviate the strongest 
from the kinematic paradigm, widely accepted for blazars, 
of outward superluminally moving jet components. BL Lacs 
appear to evolve their jet ridge lines (with respect to the 
core) less than the other source classes, hence indicating a 
slower apparent flow in their jets. On the other hand, they 
show significantly apparent wider jets than FSRQs. We have 
argued that the combination of a helically structured jet (as 
a universal property of all AGN jets) and projection effects 
would produce the apparent kinematic properties observed 
here. Viewing a source at a small angle to its jet axis pos- 
sibly allows us to uncover this peculiar kinematic behav- 
ior that would otherwise be inaccessible at larger viewing 
angles. This would then imply that BL Lac objects are pref- 
erentially seen at smaller viewing angles than FSRQs. Com- 
bined with the significant number of (non-BL Lac) CJF 
sources showing evidence that support this different kine- 
matic scheme imply a rather universal effect, rather than 
something unique to BL Lacs. 

Finally, it should be noted that the above results under- 
line the fact that the notion of linear, ballistic trajectories 
for AGN jet components usually employed until recently is 
a very crude approximation and, more often than not, devi- 
ates grossly from the reality. It is of great interest to uncover 
the physical process that leads to the properties of AGN jet 
ridge lines as studied in this paper. 
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